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Abstract 

Simple and analytically tractable expressions for functional determinants are 
known to exist for many cases of interest. We extend the range of situations for 
which these hold to cover systems of self-adjoint operators of the Sturm-Liouville 
type with arbitrary linear boundary conditions. The results hold whether or 
not the operators have negative eigenvalues. The physically important case of 
functional determinants of operators with a zero mode, but where that mode has 
been extracted, is studied in detail for the same range of situations as when no 
zero mode exists. The method of proof uses the properties of generalised zeta- 
functions. The general form of the final results are the same for the entire range 
of problems considered. 



1 Introduction 



This paper is concerned with the rather elegant, and surprisingly simple, expressions 
that exist for the functional determinants of certain types of differential operators. In 
an earlier paper P], we introduced a new method for deriving these expressions for 
operators of a relatively simple kind, which only used elementary ideas from complex 
analysis and the theory of differential equations. Here we extend the class of problems 
which may be analysed using this technique. Although the discussion necessarily be- 
comes more technical, the essential points remain the same, and we are able to derive 
the desired results without the need for any very sophisticated machinery. 

The derivation of formulae of this kind is a topic which has been investigated 
by numerous authors in the past. In our earlier paper pp, we gave a brief history 
of the subject. Essentially, most of the early results were obtained by theoretical 
physicists who were typically interested in the expressions obtained when carrying out 
Gaussian functional integrals 0-11! • These results were then extended and elaborated 
in a number of ways [Sj-jn]. However, many of these latter treatments were quite 
abstract, and also did not deal with the case where the operator has a zero eigenvalue 
(a "zero mode"). This situation is quite commonly encountered in real problems in 
theoretical physics, since in many cases a continuous symmetry in the problem is 
broken, and a zero mode is generated by Goldstone's theorem [201 • Although there 
has been some work carried out to determine the form of functional determinants 
with zero modes excluded I^Il-ESl, the methods that were used involved the use of a 
regular isat ion procedure which could have produced results which were not independent 
of the scheme adopted. 

These were the motivations for our approach described in Ref. [T]. The method 
used a generalised zeta- function [211-12^1 to calculate the functional determinants, but 
the analysis involved only methods which are familiar to theoretical physicists. It 
also covered the physically interesting situation where operators had zero modes which 
were excluded from the evaluation of the functional determinants. The method was de- 
scribed for simple operators of the type —(P/ dx'^+R{x), but for general linear boundary 
conditions. In the present paper we extend this treatment in several ways. Firstly, we 
derive the results for the general Sturm-Liouville operator —d/dx {P{x)d/dx) + R{x). 
Secondly, we allow for the fact that operators will, in general, have negative eigenval- 
ues. Thirdly, we generalise the entire formalism to systems of second-order operators. 
In all cases we derive the results for functional determinants of operators which do not 
have a zero mode, and for those which do, but where it has been extracted. 

The outline of the paper is as follows. In Section |21 we discuss the formalism for the 
general Sturm-Liouville operator, modifying our previous treatment to cover the case 
of arbitrary P{x) > and operators with negative eigenvalues. We restrict ourselves to 
operators with no zero mode; this case is discussed in Section |21 In Section 4 it is shown 
how the results of these two sections carry over to systems of r > 1 degrees of freedom. 
We conclude in Section 5 with a summary of the results of the paper and suggestions 
for future work. There are three appendices. In Appendix El we discuss the conditions 
which have to be imposed so that the operator is self-adjoint and give details of some 
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technical calculations that are required in the development of the theory. In Appendix 
IHlsome results on the asymptotic form of solutions of differential equations, which are 
used in the main text, are derived. In Appendix O some of the more technical aspects 
of dealing with zero modes in systems of differential equations are presented. 



2 One-component 

In this section we will describe our approach in the context of operators of the form 

on the interval / = [0,1]. The structure displayed in (^J is the most general that 
is possible for a self-adjoint second order differential operator of the Sturm-Liouville 
type. The functions Pj{x) and Rj{x) are assumed to be continuous on the interval 
I. In addition, we assume that the metric, Pj{x), is positive throughout the interval 
under consideration. The index j takes on only two values: the operator Li is the real 
focus of interest, but in order to control divergences in det Li, we actually consider the 
ratio det Li/ det L2, where L2 is appropriately chosen. Typically L2 will be taken to be 
"simple" in a sense that it can act as a reference with which det Li can be compared. 
The eigenproblem corresponding to (Q) is 

LjUj^Vx(^) = \u^^^{x). (2) 

Note the symmetry Uj ^ = u- It is convenient to go over to a first order formalism 
and in order to have the most natural formulation we define a new function Vj ^{x) = 
Pj{x)u'. ^{x). Then from we have that 



d / 




We will adopt the notation 



-.A-)= 5 ^.A=| ^ ^^-^ ), (4) 




in which case (jS]) may be written as 



dx 



Dj,\ix)n.^^{x) . (5) 



It is useful at this stage to introduce two unique, independent solutions of the differen- 
tial equation The solutions are made unique by specifying the "initial conditions", 
that is, the value of the solutions and their derivatives at a; = 0. Denoting these 
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two solutions by u^p^{x) and u^p^{x), the most general solution of niay then be 
expressed as 




2 / 




(6) 



We now define the two matrices 





which are related by 

It follows that det E- ,^{x) = Pj{x) det Hj ^{x). Since u^^lx) and m^.^^(x) are inde- 
pendent solutions, detHj^ ^ 0, and therefore detE'^-^ 7^ because Pj(a;) > Vx. 
Furthermore, because det E- ^[x) is the Wronski determinant for the differential op- 
erator (0)), we see that deti?^ y^(x) is independent of x. A convenient choice for the 
set of initial conditions is £^^•y]^(0) = I2 (here, and throughout the paper, is the 
m X m unit matrix). From this it follows that det ^{x) = 1 Va; G [0, 1]. Also with 

this choice for the initial conditions on u^^lx) and m^.^^(x), it follows from © that 
a = u- ^(0) and [3 = v- ^(0), that is, 

^j,A^) = ^j^Vxi^Hlxi'^) + ^,,v^(0)«i.'^(^) , (9) 
or in terms of the first order formalism. 



(10) 



So far no mention has been made of the boundary conditions on (jSj). These take 
the form of two conditions on the set ^Uj ^(0), m^. ^{0),Uj ^(1), m^. ^(1)|. These can 
be converted into conditions on u^- ^ at the boundaries, and for the case of linear 
boundary conditions 

V',,/a(0)/ V',,/a(1)/ 1,0/ ^ ' 

where M and are 2x2 matrices whose entries characterise the nature of the boundary 
conditions. Using (fTUI) these boundary conditions may be written as 



M + ATE^. ^(1) 




:i2) 
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and so the condition on A for eigenf unctions to exist is 



det 



M + iVE^.^(l)] =0. (13) 



The equations ()11|) are the most general linear boundary conditions. They fall 
naturally into two classes: 

(i) det M = 0, det = 0. In this case we can show that the matrices M and N may 
be chosen to be of the form, 



M= „ „ ; iv= " " , (14) 





that is, the boundary conditions are of the Robin type: Au- ^(0) + ^(0) = 
and Cw^. ^(1) + Dv^^^^{l) = 0. 

First, let us prove that det M = and det = implies that Mu^yi^(O) = 




and A^u^. ^(1) = 0. To see this, define u'. ^(x) = Q ^u^. ^(x) and multiply 
(jllj) by P, where P and Q are arbitrary non-singular matrices. Then defining 
M' = PMQ and A^' = PNQ, we obtain the same boundary conditions but in 
the primed system. However, since M has rank 1, we may choose P and Q in 
such a way that M' = diag(l,0) or diag(0, 1). Furthermore, since A^' has zero 
determinant it must have one of the following four forms: 

= [ka kb)^ = 

. ( a ka \ . 

= [b kb)^ 

Writing out the boundary conditions explicitly when M' = diag(l, 0) and A^' has 
each of these four forms, we find that for all cases where there are two independent 
conditions, u'.^{0) = 0, that is, M'W.^{0) = 0. Similarly, if M' = diag(0, 1), 
we find that for all four possible forms of A^', all valid boundary conditions lead 
to v[ ^{0) = 0, which for this choice of M' once again gives M'u^. ^(0) = 0. 
Returning to the unprimed system this implies that Muj ^(0) = and so from 
(fTT|) . A^u^.^(l) = 0, as required. 

We may now use the fact that M and A^ separately must have one of the forms 
Ai, . . . , y44. In each case there is only one independent relation of the Robin type. 
This may be written in the language of M and A^ matrices by adopting the forms 

m- 

(ii) det A^ 7^ 0. In Appendix B we show that this implies that det M ^ 0. Then, 
from dlH), 
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Since N is not null, neither is M~^N, and so either Uj ^{0) or Vj^{0) depend 
on the boundary conditions at x = 1. Boundary conditions such as these are 
called two-point boundary conditions, or non-separated boundary conditions, in 
contrast to the one-point or separated boundary conditions described by ()14|1 . 

After this short review of the background, we are now in a position to describe our 
method for obtaining the basic formula for detLi/ detL2. The starting point is the 
observation from (fT!?jl that the function det[M + NEjyT^{l)] has zeros at values A which 
are eigenvalues of Lj, as given by ((21) • An alternative statement is that the logarithmic 
derivative of det[M -|- NEj ^(1)] has a simple pole with unit residue at these values of 
A. This allows us to define the zeta function of Lj by 

^^^^'^ = 2^ //^^^'ZX^'''^"* K + ^^,,v/a(1)] , (15) 

where the contour 7 is counterclockwise and encloses all eigenvalues as shown in Figure 
1. 



cut for A 





A-plane 


— e-^ 


/ 











Figure 1: Contour 7 in the complex plane. 



As given, the representation is valid for > 1/2. In this section we assume that 
there are no zero modes, but note that we allow for negative eigenvalues. In order to 
avoid the negative eigenvalues lying on the cut of the complex square root, we define 
the branch cut to be at an angle 6 to the positive real axis. For most applications it 
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is the ratio of determinants of two operators that naturally occurs. This is found by 
analysing 



1 



27rz 



d , det 



/ d\\-'—\n 



dX det 



(16) 



The first idea is to deform the contour such that it encloses the branch cut of A"**. In 
order to see in which range of s-values this is possible, let us consider the large-S^-y/A 
behaviour of the integrand. 

As shown in Appendix B, we have for Pi{x) = P2{x) as $5-\/A — > ±00 the behaviour 



d det 
dX det 



M + NE,^^,{1] 



M + NE,^^,{1] 



O 



A3/2 



(17) 



So for — l/2<3?s<lwe can shift the contour such as to enclose the cut and ultimately 
we can shrink it to the cut. Taking due care of the definition of the complex root near 
the cut, we find for the upper part 



-ise 



2ni 



d det 

dX A-^4Tln 

dX det 



whereas for the lower part we have 



27ri 



det 


M + NE^^^.. 




det 







Using the symmetry ^^.^,6/2^(1) = Ej^^.n^.e/2)^{Lj, 
yield 



1), these contributions add up to 



00 



X 



dX A" 



In 



det 


M + NE^^^. 


»/2x/a(1) 


det 


M + NE^^^. 


9/2Va(1) 



sm vrs 



TT 



dX X' 



, , det 
dX det 



M + NE^ ^« 



M + NE. 
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For 9 = 71 and Pj{x) = 1 this reduces to the result of our previous paper p. This type 
of result is now perfectly suited for the evaluation of the determinant quotient. The 
prefactor disappears at s = and so 



CL,(0)-C(0) = -ln 



det [M + NEi^o{l)] 
det [M + NE2,o{l)y 



(19) 
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In particular, we note that the answer does not depend on the angle 9. Simplifying 
notation we define 




We will refer to yj{x) and yj (x) as homogeneous solutions since they all satisfy the 
equation Ljyj{x) = 0. Then the result is from (jl9p 

detLi ^ det [M + NYijl)] 
detLa ~ det [M + A^F2(1)]' 

This is formally identical to the result we obtained when Pj{x) = 1 and all of the 
eigenvalues were positive pp. This shows that this simple result is obtained even with 
the added complications of non-trivial metrics and negative eigenvalues, as long as L2 
is chosen so that P2{x) = Pi{x). 



3 Determinants with zero modes extracted 

In this section we discuss the evaluation of determinants of operators which have a 
zero eigenvalue and where this eigenvalue has been extracted in the definition of the 
determinant. We shall indicate this exclusion with a prime: thus det'L will denote 
the determinant of the operator L with the zero mode extracted. Clearly the method 
used in the last section to derive the formula for the ratio of determinants runs into 
difficulty when evaluating such determinants. Even if the contour 7 is chosen to only 
surround the non-zero values of A, the deformation of this contour to the branch cut 
will encounter the pole at the origin. Rather than dealing directly with this extra 
singularity, we can look for a function which vanishes at all the non-zero values of A, 
but not at the zero eigenvalue. This function can then be used as the basis of the 
definition of a (modified) zeta-function, from which det' L can be calculated. As we 
will show in this section, the quantity ^ = ( — 1/A) det(M + NE^ v^(l)) has these 
properties: it clearly vanishes at all the required non-zero values of A by (|T3|) and we 
will show that fi^ 7^ 0. Since we will assume that the "normalising" operator L2 has 
no zero modes, j will be set equal to 1 throughout any discussion involving zero modes. 

The first step in the derivation is relevant even if there is no zero mode. It consists 
of demanding that the solution Uj ^{x) = au^p^{x) + (3u^J\rj{x) satisfies one of the 
boundary conditions. We may choose either one of the conditions to be satisfied, but 
in general it will fix the functional form of u- ^{x) by determining the ratio of a to 13. 
In the special cases of Dirichlet and Neumann boundary conditions it will result in a 
and (3, respectively, being set equal to zero. Although the normalisation of u- ^{x) is 
obviously left undetermined, we shall now show that a suitable choice of normalisation 
results in a significant simplification of the analysis. 

To see this let us write out the boundary conditions (fTT|) in full: 

"^ll^i.VxW +"^12^j,VA(0) +^ll^i,VA(l) +^12^i,VA(l) = (22) 

m2iU^^^{S))+rn22V.^^{Q)+ri2iU.^^{l) + ri22V^^^{l) = 0. (23) 
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Now consider the explicit form of the matrix M + NE- ^(1): 

^ mil + niiu^^il) + ni2V^^''^{l) mia + niiu\^{l) + ni2vf^{l) ^ 

\ mai + n2iu^^{l) + n22t;^^^^^(l) ^22 + n2iuj.^^(l) + n22t;^^y^(l) / 



(24) 



where we have used the definition of given in ((7j). If we add (3 /a times column 

2 to column 1 of ()24|1 we get a second matrix whose first column is just a^^ times the 
boundary conditions given in ()22j) and (remembering that E- ^(0) = 12). So sup- 
pose we ask that only the boundary condition ()22p is satisfied. Since the determinants 
of both these matrices are equal, it follows that 

det(M + iVE^. ^(1)) = a-^ {m2i«^. ^(0) + m22T;,- ^(0) + n^m^^^il) + n22^;,,^(l)} 

X (-1) (mi2 + nii«5.'^(l) + n,,vf^^{l)) . 

Therefore if we make the choice 

a = - (mi2 + nnuf^^il) + ni2t;^^'^(l)) , (25) 

then 

det(M + NE.^^,{1)) = m2iM^. ^(0) + m22i;,. ^(0) + riam^- ^(1) + n22V^^^,{l) . (26) 

Similarly if we add a/ (3 of column 1 to column 2 of we get another matrix whose 
second column is just times the boundary conditions given in and (j^ . Again 
choosing the first boundary condition to be satisfied, and also now asking that ^lE^ 
holds, then we determine (3 to be given by 

13 = (mil + n^,uf^^il) + ni2.;^(;^(l)) • (27) 
So in summary, we have shown that if we take a solution of (0) of the form 

'^j,V\(^) = - {1^12 + nnuf^^il) + ni2t;^^^^(l)) ufl^ix) 

+ (mii+nii«(;^(l)+ni24'^(l))«('yx), (28) 



then 



''i":::;(o)j+n4'::.(i)i^U(M.A'^..(i))j' 

That is, if ^{x) is chosen to satisfy only one boundary condition (in this case the 
first), and its normahsation is chosen appropriately, then det(M + NEj ^{1)) will be 
directly proportional to the expression on the left hand side of the boundary condition 
(I23|) ■ with a constant of proportionality equal to unity. As far as we are concerned in 
this paper, the relation (j29|) has two important consequences: 
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(i) If there is no zero mode in the problem, we can simply take the limit A — > in the 
above formulae and get a simplified, and more explicit, expression for the result 
(1211) . To do this we choose the particular solution of the homogeneous equation 
LjUj^x) = to be 

y,{x)= - {mu + nnyf{l)+nuPAl)yf'{l))yf{x) 

+ (m^,+n,,yf\l)+n^2P,{l)yf^'{l))yf\x). (30) 



In terms of this particular solution, the A ^ limit of ()26|1 may be used to write 
the result (pT|l as 

detLi ^ m2ii/i(0) + m22Pi(0)i/l(0) + n2iyi{l) + n22Pi{l)y[{l) 
det La m2ii/2(0) + ^7122^2(0)1/^(0) + ^211/2(1) + ^22^(1)1/2(1) ' 

(ii) If there is a zero mode, instead of taking the limit A ^ 0, we use (12 (ij) as the source 
of the two relationships we need to show, namely that det(M + NE-^ ^/x{^)) ~ ^ 
for small |A|, and in particular that ^ defined earlier, satisfies /i^o 7^ 0. We 
now discuss in more detail how this is carried out. 

Let us begin by defining the Hilbert space product of ^{x) and ui q{x) on £^(/) 

by 

(^i,ol«i,v^) = / dxuifi{xyu^^^{x) , (32) 

^ 

where * denotes complex conjugation. So multiplying by ui^o{x)* and integrating 
gives 



1 

dxuifi{x)*LiU-^^y^{x) = A (mi,o|mi,v^) 



By partial integration we get boundary terms plus LiWi o(x)*. This latter term is zero, 
so therefore 



Ml 



We can now use to solve for two members of the set Va(0), "^i v^(0), Ui ^^(1), 

^(1)| in terms of the other two and det(M + NE-^ ^(1)). An exactly analogous 

procedure is carried out on the set {mi^o(0)*, fi,o(0)*, Mi,o(1)*, "^1,0(1)*}) but in this case 
Uifi{x) satisfies both of the boundary conditions, and so (fTT|) . rather than (|^. should 
be used. This procedure is discussed in more detail in Appendix A, where it is shown 
that substituting the expressions for these four quantities into the left-hand side of 
eqn. (jHHj) shows that it is directly proportional to det(M + NE-^ v^(l))- The constant 
of proportionality (denoted by B~^) is independent of A, and only depends on the 
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nature of the boundary conditions and on the A = solution ui^o(a;) at the boundaries. 
Therefore we may write 



det(M + NE,^^{1)) = B k^^(x)t;i,o(x)* - u.^oixYv.^^.ix)]' = S A K,o|%,va) -(34) 



It should be stressed that while uifi{x) satisfies both boundary conditions, y^i^) (A 7^ 
0) satisfies only one (together with a normalisation condition), in other words, it has 
the form ()28|) . If the other boundary condition is imposed, A is restricted to take 
on values for which A is an eigenvalue, and we see from (jH^ that the orthogonality 
condition {uifi\u^ ^) = holds by virtue of (fT^ . 

The constant B is determined in Appendix A, where the related question of the 
conditions for the operator to be self-adjoint, is also discussed. The conclusions are: 

(i) If the boundary conditions are separated, so that det M = det = 0, then if the 
operator is self-adjoint, M and A^ can always be chosen to be of the form ()14p. 
with M and A^ real. In this case 

^ = ' if ^21 ^ ; B = ^ , if n22 ^ 0. (35) 

^^i,o(l) wi,o(l) 

(ii) If the boundary conditions are non-separated, so that det M 7^ 0, det A^ 7^ 0, 
then if the operator is self-adjoint, M and A^ can always be chosen so that one 
of them is real, say A^ = N^, and the other one a real matrix times a phase: 
M = Mrc'"', < a < 27r. It also follows that det Mr = det Nr. In this case 

ni2n2i - niin22 

^^11^^1,0(1) + ^12^^1,o(l) 

The function ^ mentioned earlier can now be identified. If we define 
_ det{M + NE,y-,{l)) 

/l,^/A = ^ ■ = (^1,0|Mi,Va) ' (37) 

we see that it vanishes at the required values of A, but is non-zero when A = 0. However, 
in our evaluation of the contour integral in the last section, it was also vital that the 
large |a/A| behaviour for j = 1 and j = 2 were the same, so actually we need to replace 
det(M + NE-^ Va(1)) the integrand of the contour integral by (1 — A)/^ This has 
the required properties when both A = and A 7^ 0, but in addition it behaves like 
det(M + A^i?^ Va(1)) large |V^|, also as required. So in order to derive an expression 
for det' Li/ det L2 we need to begin from 



27ri J-y dX det 



where the contour 7 encloses the point A = 1 and the values of A on the real axis which 
define the eigenvalues. 
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It is understood that the zero mode has been omitted from the definition of the 
zeta function. For definiteness we have assumed that the contour encloses A = 1 such 
that the term '-1' on the right hand side corrects for the contribution due to the factor 
(1 — A). Proceeding as before, now noting that e^s/^v^ = /i,e-»{''-«/2)VA' obtain 



oo 



ClAs) - ClM = -1 + e-(--^)^ / dX A-^^in 



d , (l+e^'A)/i_,..p^ 





For the derivative at s = this means 



dX det{M + NE^^^,e/2^il))' 



Using the notation of equations pO|l and (|37jl . this may be cast into the final form 

det'Li B{yi\yi) 



detL2 det [M + A^r2(l)]' 

4 Systems of differential operators 



(39) 



The extension from a single differential operator of the form to a system of differ- 
ential equations is relatively straightforward, the main problem being one of notation. 
Provided that the previous sections have been read, the discussion in this section should 
be clear, since it parallels the case of a single operator. As an additional aid to under- 
standing, we illustrate new concepts which are introduced on a specific example. Some 
of the more cumbersome formulae which are not vital to an overall understanding of 
the formalism are relegated to Appendix O 

We consider the system of differential operators 

where Rj{x) is a Hermitian r x r matrix and j = 1, 2 labels the two different determi- 
nants. We assume Pj{x) to be scalar, which is the relevant case for most applications. 
The second order problem is rewritten as a first order problem in the standard way. 



with the matrix 



3 

Rj — X • If Oj. 



and where now u- ^(x) and ^(x) are r-dimensional vectors. We define, as before, 
the fundamental matrix as 




12 



with u^^^{x), cr = 1, . . . , 2r, being again r-dimensional vectors. The boundary 

conditions read 



or, alternatively 



So the condition for the eigenvalues reads 

det(M + iVi?^.^(l)) = 0. (42) 

In the case that Pj{x) is scalar, the analysis in Appendix B goes through. The only 
change is that the heat kernel coefficients contain a trace over the internal degrees 
of freedom. With this change, the asymptotic behaviour of the relevant integrand is 
known and for Pi{x) = P2{x) we can proceed as previously. In the absence of zero 
modes we find formally the same answer as before, 

detLi _ detjM + NYijl)) 
det L2 ~ det{M + NY2{1))' 

Let us next consider the case with zero modes. In order to explain the individual 
steps of the general formalism that we are developing, we will illustrate each step using 
a specific example encountered in the study of transition rates between metastable 
states in superconducting rings |2Z1I2H]- The differential operator in this problem is 
defined on the interval [— //2, //2] and has the form 



V (1 - /i2)e-2^^- _^ + (i_2/i2) 

Boundary conditions imposed are so-called twisted boundary conditions defined through 

M = -diag(e^'^', e"^^', e^'^', e"*'''), N = J4. 

We will refer back to this example at suitable stages of our procedure. 

The starting point for the general formalism is as before, namely ()33|). If we can 
derive a relationship of the form 

det(M + NE,^^,{1)) = B[u,^^,{x)v,,oixr - «i,o(x)*t;i,^(x)]J , (44) 

where B is known, then from (|33|) we have that 

det(M + iVE^,^(l)) = B\ {u,,o\u,^^) . (45) 



= -^^2 + i?i(x). (43) 
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This is precisely as in Section and allows us to identify the function defined 
by (|H7|). which is to be used in the proof of the result. 

So, let us return to the proof of (j44|) . We will show that, if we appropriately 
normalise Uj ^{x), then by imposing all of the boundary conditions but one — so that 
A is not constrained to be an eigenvalue — we may write 



M 



^i,Va(0) 



+ N 



^i.v^(l) 








det{M + NE,y-,il)) ) 



(46) 



This equation is exactly analogous to in Section El where we imposed only one 
out of the two boundary conditions. Here there are 2r boundary conditions and we 
will impose 2r — 1 of them. 

To obtain we first write Uj ^{x) as a linear combination of the 2r fundamental 



solutions 



2r 



U„ 



2r 



X] 



X] 



(47) 



(T=l 



CT=1 



where we have dropped the j and dependence from the a. Since E- ^(0) = /2r, 



a 



(^) 



u 



if a = 1, . 



^^_,(0), if a = r + l,...,2r^ 



where Uj^^{x) is the ath entry of the vector u-^{x), with a similar notation for 
^{x). Therefore using (j4H) . but only imposing the first 2r — 1 boundary conditions. 



3 

gives 



(M + iVE^.^(l)) 



a 



(2) 



/ \ 




V * J 



(4J 



First, suppose that det(M + iVE^. ^(1)) ^ 0. Then, multiplying (gHl) by (M + 
iVE^.^(l))-i yields 



adj(M + iVij;^.^(l)) 
det(M + iVE^.^(l)) 



/ \ 




V * / 
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where adj(M + NEjyy^{l)) is the adjoint of the matrix M + NEj^^{l). We see that 
the choice of det(M + NEj^{l)) for * is natural, since in this case the expansion 
coefficients have the simple form 



a^'^) = adj(M + iVE,^(l)).2r-. 



(49) 



In the r = 1 case this simply leads to the results ^ and If det{M+NEj^^{l)) = 
0, then by (jl^ Uj^{x) is an eigenfunction which satisfies the boundary conditions, 
and so (jlFIjl also holds. 

Altogether there are 4r boundary data, r data coming from each of u^^{0), 
"^iVaIO)' '^l^/x(^)■ ^'i- allows us to express 2r of the boundary 

data in terms of the other 2r data, which we call the complementary ones. Suppose 
that 6 is a vector consisting of the 2r boundary data that we wish to express by the 
complementary ones, collected in be- Expressing pUj) in terms of these values gives 



Zb 



Zcbc 








\ 



1^ det (m + NE.^^^il)) I 



(50) 



where Z and Z^ are (2r x 2r) matrices built from the various components of M and A^. 
To state b,bc, Z and Zc explicitly, we need to introduce several indices which refer to the 
ways in which the boundary data are re-distributed within each of the four boundary 
data groups. Let i, j, k, I be indices all of which can take on values from to r and such 
that i + j + k + I = 2r. Let {ai, a,.} and {ci, c^} be arbitrary permutations of 
the numbers {1, ...,r}, and also let ...,br} and {di, ...,dr} be permutations of the 
numbers {r + 1, ...,2r}. These index groups are such that m^^ acts on boundary data 
in ^(0), rrih. acts in ^^(0), ric^ acts in ^^(1), and finally n^^ acts in ^^(1)- 
The general form of b,bc,Z and Zc are discussed in Appendix O from which it is clear 
that b can be expressed through be only if the matrix Z is invertible. The choice of 
the data b has to guarantee this is indeed the case. That this is always possible follows 
from the fact that M and A^ define boundary conditions such that has a unique 
solution for A an eigenvalue. If a suitable choice of b were not possible, the boundary 
value problem would not have a unique solution. 

For the example described by the most natural choice for b, be, is 



/ «l,v/A,l(V2) \ 

«i,Va,2(V2) 

^l,v/A,l(V2) 

V ^i,v^,2(V2) J 



( «i,vO:,i(-V2) \ 

«1,VA,2(-V2) 
^l,v/A,l(-V2) 
V ^l,v^,2(-V2) j 



(51) 



so that Z = N (= Ia) and Zr 



M 



-diag(e*'^', e «A''^e*^',e *^')). This guarantees Z 



is invertible and b can be expressed through be- 



15 



Alternatively one could, for instance, choose 



^(alt) 



/ «i,VA,i(-V2) \ 
«i,Va,2(-V2) 

^l,v/A,l(-V2) 

V ^i,v^,2(-V2) j 



^(alt) 



/ «l,VA,l(V2) \ 

^i,Va,2(V2) 
^i,v^,i(V2) 

V ^1,V^,2(V2) ) 



(52) 



In this case 



^(alt) ^ ^(alt) ^ ^_ 



Again, Z is invertible and 6 can be expressed through he- Clearly, there are many other 
choices of 6, he and the associated Z^. 

Going back to the general formalism, given a suitable particular choice of Z, this 
allows us to express the 2r data h by the complementary 2r data he- The explicit 
expression is given by equation ()(^3jl in Appendix O The entries of h can now be 
substituted into the left-hand side of (jHHj) and the terms collected together. As discussed 
in Appendix O this leads to (IHj) with 



^ ^ — X! ^(i+i+a)(2r)'^l,0,CQ(l)* ^ ^ {i+j+k+a)(2r)'^^fi4a-riX)* 



a=l 



^a(2r)^l,0,a,(0)* + Yl ^(i+a)(2r)"l,0,6c-r(0) 



(53) 



a=l 



where Z^^}^ refers to the (/37)-component of Z~^. 

To illustrate the use of this result let us apply it again to the example 
For the choice ()5H) we have i = 0,j = 0,k = 2,l = 2 and we obtain 



KA^>i,oi^r - ^ifiixryi,M^)]-i/2 = - det(M + ArE,_^(//2))ni,o,2(V2)*. 
For the choice ()52j) we have i = 2,j = 2,k = 0,l = and we obtain 

K,VA(^)^i,o(a^)* -«i,o(a:)*t;i,^(x)]L/?/2 = -det(M + ArE^^^(//2))e^^V,o,2(-V2)^ 



Comparing with (j4^ we see that for the choice (j5T| B^^ = — mi^o,2(V2)* and for the 
choice (IS2I) B^^ = — e*^'Mi_o,2(— V2)*- Taking into account the boundary conditions for 
the zero mode ui,o{x), the two answers are seen to agree. Furthermore, this answer 
agrees with the result calculated in pSj . 

Returning to the general case, we see that we have proved the result (jlSj) with B 
given by (j^Hj) . The proof now proceeds as in Section El and we once again find the result 
(jSni)- The function yi{x) in this result is the zero mode, and satisfies the boundary 
conditions, but it has to be appropriately normalised: 



2r 



X) 



(54) 



<7=1 



where the y['^\x) are the 2r fundamental solutions chosen to satisfy 1^1(0) = I21 



2r- 
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5 Conclusion 



The two main results of this paper are ()2H1 and They give expressions for the 

ratio of functional determinants in terms of the nature of the boundary conditions and 
the solution of the homogeneous equations formed from the operators in question. The 
first result holds if the equations have no zero modes and the second holds if such an 
eigenvalue exists, but has been excluded from the evaluation of the functional deter- 
minant. These results agree with those obtained in a previous paper pLj, but now the 
range of operators for which they are valid have been considerably extended to: those 
which have a metric Pj{x) ^ 1, those with negative eigenvalues and systems of opera- 
tors. Although the results are simple to state, a slightly more thorough appreciation 
of the method is required in order to apply them to a particular case. For instance 
the solution yi{x), which appears in the results, is the solution of the homogeneous 
equation satisfying the boundary conditions. This solution is only defined up to a con- 
stant, but a particular choice for this constant has to be made if the simpler form of 
(PT|) — given by — or the zero-mode result is to be used. An explicit form 
for yi{x) is given by ()54|). The choice of normalisation originates from requiring that 
the right-hand side of ()29|) or pUj) is the required determinant, but with a constant of 
proportionality which is equal to 1. 

Once the suitably normalised solution yi{x) has been obtained, the rest of the cal- 
culation is straightforward. The result only depends on this function — and on none of 
the other eigenfunctions — and on the matrices M and which define the boundary 
conditions of the problem under consideration. In addition, in some applications, the 
norm {yi\yi) will cancel out with the Jacobian of the transformation to collective coor- 
dinates, and therefore will not be required. In this case, however, it will be necessary 
to check that the zero-mode has the same normalisation as has been adopted in the 
derivation of (jH^ . So for this case only the properties of yi{x) at the boundaries would 
be required. If an analytic expression for yi{x) cannot be obtained, there should be 
little difficulty in obtaining numerical values for the boundary data on this function. 
While the proof of the results which we have obtained are easily accessible, they need 
not be understood in order to apply the results to a particular problem. 

We believe that the results presented here cover a wide range of problems where 
they are likely to prove useful. There are still a number of possible extensions that 
are open to investigation. Examples include operators with derivatives higher than the 
second, single determinants rather than ratios of determinants, and operators in more 
than one dimension. We hope that, in addition to the concrete results which we have 
obtained, this paper will serve to stimulate work on these and related problems. 

Acknowledgements: The research of K. Kirsten was partially supported by the 
Max Planck Institute for Mathematics in the Sciences (Leipzig, Germany) and the 
Baylor University Summer Sabbatical Program. 
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A Self-adjoint condition and related questions 

In this appendix we will consider two technical points encountered in Section IHl They 



are 



1. Condition for problem to be self-adjoint. 

The boundary conditions considered in section 121 are given by ()22|1 and ((2 



m2iu.y^{0) + m22Vjy^{0) + n2iu.yj^{l) +n22V-y^{l) = 0. (Al) 



Suppose that U^^l^i^) 

and f/j^^( x) are any two functions (which are not, in 
general, solutions of (j21)) which satisfy these boundary conditions. The condition 
for the problem to be self-adjoint is that 

\uZd-)V;^{xr - U^-rv;^ix)\ 1 = 0, (A2) 

where, as in the main text, Vj ^{x) = Pj{x)U'. ^{x). We wish to solve for any 

two members of the set {L^j^^(O), K^^^(O), f/j^^(l), K^^^(l)} in terms of the 
other two by using the two boundary conditions (and similarly for the second 
solution II). Substituting these four functions into ()A2|) in terms of the other 
four will give us the conditions that need to be imposed on the matrices M and 
for the problem to be self-adjoint. 

2. Proof of the first equality in eqn. \34^ . 

In this case the A = solution Uifl{x) will satisfy the boundary conditions ()A1|) . 
but the A 7^ solution will only satisfy one boundary condition and a normali- 
sation condition, that is (see eqn. (j^ ): 

"^ii«i,Va(0) + "^12^^1,^(0) +^i2^i,Va(1) = (A3) 
'^2i%,VA(0) + '^22t^i,v^(0)+n2in,y^(l)+n22t;i,^(l) = det{M + N E ^ ^^{l)) 

As discussed in Section El we wish to solve for any two members of the set 
{■"i Va(0)? "^1 ^/a(0)) ""i ^/a(1)) "^1 v^(l)} ™- terms of the other two and det(M + 
NE^ v^(l)) by using the conditions ()A3|) . Similarly, we wish to solve for any two 
members of the set {mi_o(0)*, fi^o(O)*, Mi,o(1)*, "^1,0(1)*}; in terms of the other two, 
this time using ()A1|1 . Substituting these four functions into the left-hand side of 
in terms of the other four, will enable us to show that 

«i,VA(a;)^i,o(a;)* - wi,o(x)*t;i,v^(x)]^ (X det(M + NE^^^{1)) , (A4) 

where the constant of proportionality is independent of A, and only depends on 
the nature of the boundary conditions and on the A = solution ui o(a;) at the 
boundaries. 
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It is clear that these two questions are related. In fact, the proof of the first point is a 
special case of the proof of the second; we simply need to set det{M + NE^ ^/a(1)) equal 
to zero everywhere. We will prove the result first in the case of separated boundary 
conditions and then for non-separated ones. 



For separated boundary conditions we have det 
and may be chosen to have the form (fT^ . 



and det M = 0, and M 



If mi2 7^ and n22 7^ 0, then ()A3jl may be written in the form 

mil 



^i,vO:(0) 



-u 



i,Va(0) 



mi2 



n2i 



-u 



1,VA 



11 



n22 n22 

Equivalent results hold when A = if the determinant is set equal to zero. 
Eliminating the f 's in terms of the u's yields 



n22 



det(M + ArE^^^(l)) 



+ «i,o(1)*Wi,Va(1) 



n2i 
n22 



n 



21 



n. 



22 



«i,o(0)*w,^^(0) 



mil 
mi2 



m 



11 



m 



12 



If we first of all assume that all of the boundary conditions are satisfied, then the 
determinant is not present and we see that the general condition for the operator 
to be self adjoint is that the ratios mii/mi2 and n2i/n22 be real. Since we always 
have the freedom to multiply the first line of ^\\\ by an arbitrary complex number 
and the second line by another arbitrary complex number, we can always choose 
mi2 and n22 to be real, in which case we deduce that mn and n2i should also be 
real. Therefore if the operator is self-adjoint, the matrices M and N can always 
be chosen to be real. 

If mi2 = 0,n22 7^ 0, then ^i ^(0) = and mi,o(0) = 0. The above expression 
then tells us only that n2i/n22 must be real if the operator is to be self-adjoint. 
But now mil and n22 may be chosen to be real, and we once again find that M 
and may be taken to be real. The remaining cases where n22 = may be 
treated in the same way. 

In summary, when det M = det A^ = 0, if the operator is self-adjoint then M and 
A^ may always be chosen to be real, and 

1 





"i.o(l)* 

"22 



det{M + NE^y^il)), ifn22^0 
detiM + NE.^^il)), if nsi ^0 . 



(A5) 



Through the boundary condition for the zero mode, these two forms are clearly 
equivalent if both ^21 and n22 are non-zero. 
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(ii) Suppose that det 7^ 0. Then multiplying by ^ and taking j = I gives 

/%,v^(l)\ ^ -ni2 W 

U,vO:(l)/ detiV 1^-^21 nu ) [det{M + NE^y-,{l)) 

( du d,2 \ KA^)] 

\ d21 d22 ) [v.^^MJ ' 

where the dij are the elements of the matrix D = N^^M. Substituting for 
u^^{l) and fj^^(l) (and their A = counterparts, which do not contain the 
det{M + NE^^^{1)) term) gives 

Ui,^{x)vifi{x)* - Mi,o(a;)*t^i,^(a;)]^ 
^ _ (nn^i,o(l)- + n,2Vr,om ^^^(^ ^ 

+ Mi,o(0)*Wiyi^(0) {1 + (^12^21 - dl^d22} + t;i,o(0)*t;iyx(0) {^126^22 " ^^12^22} 

- ^^i,o(0)*Miyx(0) {1 + d*i2d2i - c?ii'i22} " Wi,o(0)*Mi^^(0) {rf2ic?n - c/gic/n} . 

If all boundary conditions were satisfied, the first term on the right-hand side of 
this expression would be absent, and the general conditions for the operator to 
be self-adjoint are given by the vanishing of the brackets in the four remaining 
terms: 

diid22 ~ C^12'^21 ~ l ) '^ll'^21 ~ d*id2i 

dudl2 — C?12'^21 = 1 ; C?12C?22 = '^12'^22 • (A6) 

Examination of the conditions ()A6|) shows that they may be written in the alter- 
native form 

dij = rije^°, riir22 - ri2r2i = 1, < a < 27r, r^j G IR. (A7) 

In other words, D = i?e*", where i? is a 2 x 2 real matrix with entries rij and 
deti? = 1. Therefore if we multiply (jllj) by A^~^ we see that we may take 
M = i?e*", N = I2, OT equivalently if we multiply by a real non-singular matrix 
Nr, M = Mrc'", N = Nr where Mr = NrR. Note that det Mr = det Nr. 
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B Asymptotic behaviour of solutions at endpoints 



In this appendix we are going to analyse the behaviour of (d/dX) In det [M + A^ii^^- 

for Q^\/A — ±00 as it is needed in equation (fTHjl . We will omit the index j and consider 
the general Sturm- Liouville problem 

with the boundary conditions (lllj) imposed. The zeta function associated with this 
problem is then given by equation (fT^. 




^^^'^ = 2^/"^^ A-^lndet[M + iVi?^(l)]. 



The meromorphic structure of the zeta function is determined by the large-9V^ be- 
haviour of the integrand. The results in |29j suggest that as ^^/X — > ±00 the asymp- 
totic expansion has the general form 

d °° 

-lndet[M + iVE^(l)] ~ Y.i^r^'^n-i. (Bl) 

"'^ n=l 

Note that exponentially small terms have been dropped. 

We will now show that the coefficients An-i are related to the associated heat kernel 
coefficients and use this correspondence to prove that the first two coefficients do not 
depend on R{x). As a consequence we can trivially conclude the behaviour (|T7jl . 

We start summarising some well known facts about the heat kernel coefficients and 
their relationship to the zeta function [SOIES!- The heat trace is defined as 



where A; are the eigenvalues of the operator under consideration. As t — > 0, this sum 
clearly diverges, since we are summing over infinitely many eigenvalues. The behaviour 
as t ^ may be extracted from a classical theorem of Weyl [SI] , which, in the present 
context, states that for a second order elliptic differential operator the eigenvalues 
behave asymptotically for / ^ 00 as 



7ll 



' /o 



With the help of a resummation, 



Z=— 00 l=—oo 
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it is seen that this imphes K{t) = 0(t ^/^). In more detail one can show the asymptotic 
t — s> behaviour 

oo 

i^W~E«/^"'^^'> (B2) 

j=0 

where exponentially small terms as t — > have been neglected. Here, ai are the so- 
called heat kernel coefficients. They depend on P{x), R{x), and on the boundary 
conditions imposed. We have, for example, 

Oo = (47r)~^/M dx , ai = c{M,N), 

Ip(x) 



where the constant c(M, A^), as indicated, depends on the boundary condition imposed. 
The next coefficient 02 involves the dependence on R{x). As this is of no relevance for 
us, we do not display higher coefficients. 

The heat kernel coefficients determine the residues and certain function values of 
the zeta function. To show how the relationship is derived we assume that no zero 
modes are present; otherwise, in the following calculations we have to exclude them 
explicitly. 

First by definition 

1 /"OO 1 rca 

f^o r(s);^Jo T{s)Jo 

valid for > 1/2. As is clear, the meromorphic structure of Cl{s) is related to the 
t ^ behaviour of K{t). Thus the poles of CL(s)r(s) are determined by the integrals 



0-l)/2 

j=0 



In detail we have 



Res Cl{z) = ~~Y~' ^ ^ 

Cd-q) = {-iyq\ai+2g for g G INq. (B3) 

The asymptotic expansion ()B1|) determines the above properties of the zeta function 
and thus relates An and a„. Proceeding as before, see (fTHj) . we shrink the contour to 
the branch cut at the angle 6. For the case without zero modes, as A ^ we have the 
behaviour A~'^ and as A ^ 00 we have A~'^~^/^. The A behaviour imposes 3?s < 1, 
whereas the X —>■ 00 behaviour imposes 3fts > 1/2. This shows, that the representation, 
as given, is valid for 1/2 < 3fJs < 1. It also shows, that the residues and function 
values, ()B3|) . which all lie to the left of 3fts > 1/2, are solely determined by the large- A 
behaviour. Keeping only the relevant terms to reproduce ()B3|1 . we continue 



e 



is(TT-e) sin(7rg) ^ ie(i-n/2) An-\ 



ti s-l + n/2' ^ ^ 
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which can be analysed easily in the whole complex plane. We see that for n odd, 



which shows 



whereas for n even, 



and so 



Res Cl ( 1 - - j = -Arr-l, 



= — f(r!^, (B5) 



/l„_i = (-l)"/2-i(^/2-l)!a„_i. (B6) 



In particular. 



Aq = — = -- /^rfx ^^=, Ai = ai = c(M,N), 

r(i/2) 2 Jo J]^y ^ ' ^' 



and (fT7|) follows. 

Note, that in ()B4jl we have e*^/^-\/A with positive imaginary part. A negative imag- 
inary part, such as in — e*^/^-\/A, changes the sign of A„_i for n odd. 

If there are zero modes, say r in number, the equation for ai changes slightly. First, 
given we exclude the zero mode from the definition of the zeta function, we have now 



1 °° 1 



This shows, that ()B3|) remains unchanged apart from 

a(0) = ai-r. 
Given there are r zero modes, for |A| <^ 1, we have 

±lndet[M + NE^,{l)]^j + ... 



Repeating the discussion below ()B3|) . we see that this time the A ^ behaviour 
imposes < 0, which contradicts the condition 3?s > 1/2 from |A| oo. Therefor, 
we cannot shrink the contour to the cut, but instead use the contour given in Figure 
121 consisting of a small circle 71 of radius e, and of 72 being the part of 7 shrunk to the 
cut. 

Along the contour 72 we can proceed as previously and obtain 



00 



a,,.(3) = e-(--^)^^^ J dX A-^^lndet[M + ArE^,«/.^(l)]. (B7) 
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A-plane 



71 



Figure 2: Contour 7 in the complex plane. 



For the contributions along the circle 71 we obtain 



= -e"<-)!^r^. (B8) 

TT S 



The contribution of ()B7|1 to the quantities in ()B3|1 are evaluated precisely as before. In 
addition, as s — > 0, ()B8|1 produces 

a,,,(0) = -r. 

As a result, the asymptotic behaviour is determined again through equations ()B5j) and 
(El. 



C Zero modes in systems of differential operators 

In Section 0] we introduced the vectors b and and the matrices Z and Z^. The 2r 
dimensional vectors b and 6c together contain the 4r boundary data: r data coming 
from each of u-^ ^^(0), ^^(0), ^(1) and ^^(1)- The 2r x 2r matrices Z and 
together contain all the elements of the matrices M and but rearranged in a way 
which corresponds to the organisation of the boundary data in b and be- The purpose 
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of this appendix is to make explicit the notation required to describe which of the 2r 
boundary data goes into b and which goes into be and which of the elements of M and 
N go into Z and which go into Zc. 

To do this we introduce indices i,j,k,l and permutations {ai, ...,ar}, {bi,...,br}, 
{ci, Cr} and {di, dr} as in SectionlH Let us recall that these index groups are such 
that nia^ acts on boundary data in u-^ v^(0), acts in vO^(O), n^^ acts in ^{^), 
and acts in ^^(1)- So if 

( «l,v/A,a,(0) \ 

"l,VA,a,(0) 
^l,VA,f,i-r(0) 

^l,VA,6,-r(0) 
^l,VA,ci(l) 



U 



l,vO:,Cfc(l) 



, be 



V ^l,v^,d,-r(l) / 



,v^,6j+i-r(0) 



^l,v^,^-r(0) 

^i,Va,c.(1) 

^l,VA,d,+i-r(l) 



then 



and 



Z 



Zr 











\ 'm{2r)ai ■ 


■■ m^2r)a. 


m(2r)hi ■■ 


• 'm(2r)bj 




■ rile. 


nidi 


nidi \ 


^(2r)ci •• 






n{2v)di J 




■ ■ mia. 




■ ■■ mib^ 








■■■ ni(^2r)br 




■ ■■ nie. 


nidi+i 


■ ■■ nid. 




■■■ 'n{2r)cr 


n(2r)di+i 


■■■ n(2r)dr 



(CI) 



(C2) 



The notation is such that if one of the indices i,j,k,l equals zero, then the corre- 
sponding entries above are simply absent. It is clear from (jSUjl that b can be expressed 
through be only if the matrix Z is invertible. If this is the case then, for a particular 
choice of Z, (j5(jp allows us to express the 2r data b in terms of the complementary 2r 
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data be as 



/ 


"l,VA,ai(0) 






' «l,v^,a.+i(0) 


\ 




/ \ 




«l,v^,a,(0) 






«l,V^,a,(0) 
















^l,VA,b,+i-r(0) 










^l,vO:,fe,-r(0) 






^l,v^,fe.-r(0) 


























"l,v^,Cfe(l) 
















^l,VA,di-r(l) 















V 


^l,VA,d,-r(l) 












1^ det(M + iVE,^^(l)) ) 



(C3) 



We may now substitute the 2r quantities on the left-hand side of ()(y3|l into the left- 
hand side of (jnSl)- As we have argued for the case r = 1, all terms that do not depend 
explicitly on the term containing det(M -|- NE^ yA(-'-)) have to cancel each other due 
to the self-adjointness of the boundary value problem. Without attempting to state 
the conditions required for M and iV (this does not do any harm simply because they 
are not needed), we keep only terms that do depend on the term det(M + NE^^^{\}j)^ 
knowing the others have to cancel. In this way we arrive at 

Kv^(^Ko(a;)*-Wi,o(x)Xv^(^)]o = det(M + iVE,y3^(l)) X 

( k I 

{ 51 '2^(i+j+a)(2r)^l,0,Ca(l)* ~ X! ^(i+j+fc+a)(2r-)'"l>0,da-r (1)* 

- E Z-l2rf,,o,aA^r + 2(7+«)(2.)«lAfe.-r(0)* [ , (C4) 
a=l a=l J 

where Z^^ refers to the (/?7)-component of Z^^ . This is of the desired form ()44|) with 
B given by (f^ . 
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